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Abstract. We consider tlic A^-body problem in spaces of constant curvature 
and study its rotopulsators, i.e. solutions for which the configuration of the 
bodies rotates and, usually, changes size during the motion. Rotopulsators fall 
naturally into five groups: positive elliptic, positive elliptic-elliptic, negative 
elliptic, negative hyperbolic, and negative elliptic-hyperbolic, depending on the 
nature and number of their rotations and on whether they occur in spaces 
of positive or negative curvature. After obtaining existence criteria for each 
type of rotopulsator, we derive their conservation laws. We further deal with 
the existence and uniqueness of some classes of rotopulsators in the 3-body 
case and prove two general results about the qualitative behaviour of these 
orbits. An interesting finding is that of a class of rotopulsators that behave like 
relative equilibria, i.e. maintain constant mutual distances during the motion, 
but cannot be generated from any single element of the underlying subgroup 
SO{2) X S0{2) of the Lie group SO{4). 



1. Introduction 

The goal of this paper is to study rotopulsators of the curved N-body problem, 
i.e. a type of solution that extends the concept of homographic orbit from the Eu- 
clidean space to spaces of constant curvature. In curved spaces, similar geometric 
figures are also congruent, so the word homographic is not suited for describing 
orbits that rotate and expand or contract during the motion. We will therefore 
introduce the concept of rotopulsating orbit, which coincides with the notion of 
homographic solution when the curvature of the space becomes zero. But before 
describing the results we prove in this paper, let us give the overall motivation 
of this research and provide a brief history of the iV-body problem in spaces of 
constant curvature. 
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2. Motivation 

The curved A^-body problem offers an opportunity to look into the nature of 
the physical space. How do we measure the shortest distance between two points: 
along a straight line (like on a fiat table), along an arc of a great circle (as between 
Paris and New York), or along a geodesic of some other manifold? Apparently, 
Carl Friedrich Gauss tried to answer this question by measuring the angles of 
triangles formed by three mountain peaks, to decide whether their sum was smaller 
or larger than vr radians. But his experiments failed because no deviation from 
71 could be detected outside the unavoidable measurement errors. It thus became 
clear that if the physical space is not flat, the deviation from zero curvature is very 
small. Some recent cosmological experiments involving the background radiation 
made physicists believe that space is Euclidean, although they have no proof so 
far, their results being at this point as inconclusive as those of Gauss. In fact, from 
the mathematical point of view, zero curvature is highly unlikely, if compared to 
a continuum of possibilities for positive or negative curvature. 

The study of the curved A^-body problem, however, may offer some new insight 
in this direction. We cannot travel to distant stars to measure angles of triangles, 
but we can observe celestial motions from Earth. If we prove the existence of orbits 
that characterize only one of the positive, zero, or negative constant curvature 
spaces, i.e. don't occur in any of the two others, and are stable too, then we can 
hope to decide the shape of the physical space by astronomical observations. 

This new approach towards understanding the shape of the universe is likely to 
succeed at the local, but not the global, level because celestial motions at large 
scales are mainly radial: galaxies and clusters of galaxies merely move away from 
each other. Interesting celestial orbits occur only in solar systems. But the study 
of the curved A^-body problem offers a new mathematical play field that can shed 
some light on the Euclidian case, by studying the bifurcations that occur when 
the curvature tends to zero, and may lead to a better understanding of several 
mathematical questions, including that of singularities in the motion of particle 
systems, [12], [I]- 

3. A BRIEF HISTORY OF THE CURVED N-BODY PROBLEM 

In the 1830s, Janos Bolyai and Nikolai Lobachevsky independently proposed a 
2-body problem in the hyperbolic space H^, ||2j, [T7]. They suggested the use of a 
force that is inversely proportional with the area of a sphere of radius r, where r 
is the distance between the bodies. Their line of thought followed that of Newton, 
who had viewed gravitation as inversely proportional to the area of an Euclidean 
sphere. But neither Bolyai nor Lobachevsky came up with an analytic expression 
for this new force. In 1870, Ernest Schering pointed out that, in H^, the area of a 
sphere of radius r is proportional to sinh^ r, so he defined a potential that involves 
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cothr, Wilhelm Killing naturaly extended this idea to the sphere with 

the help of a potential proportional to cotr, fT3]. 

In the first years of the 20th century, Heinrich Liebmann studied the Kepler 
problem (the motion of a body around a fixed centre) for the cotangent potential 
and recovered Kepler's laws in slightly modified form, [15] . Moreover, he found 
that a property earlier proved by Joseph Bertrand for the Newtonian potential, 
namely that all bounded orbits of the Kepler problem are closed, [1] , was also true 
for the cotangent potential, [16]. These results, together with the fact that both 
the Newtonian and the cotangent potential of the Kepler problem are harmonic 
functions in the 3-dimensional case, established the cotangent potential among 
the problems worth researching in celestial mechanics. More recently, results in 
this direction were obtained by several Spanish and Russian mathematicians, [3], 
[H], [22] • It is interesting to note that, unlike in Euclidean space, the curved 
Kepler problem and the curved 2-body problem are not equivalent. While the 
former is integrable, the latter is not, [22], so its study appears to be complicated 
than that of the classical case. 

Recently, we found a new setting that allowed the generalization of the 2- 
body case to any number > 2 of bodies, [TT], [6], [7], where we showed that 
the equations can be simultaneously written for positive and negative curvature. 
The idea was to use the hyperbolic sphere, i.e. Weierstrass's model of hyperbolic 
geometry, given by the upper sheet of the hyperboloid of two sheets, embedded 
in the Minkowski space. By suitable coordinate and time transformations, this 
study can be reduced to and H^, or to and in the 2-dimensional case. For 
the latter, the equations in intrinsic coordinates were also obtained, [20], [ID]- So 
far, we studied the singularities of the equations and of the solutions, the various 
classes of relative equilibria, and some rotopulsating orbits in the 2-dimensional 
case, [1], [5], [6], [7], [9], [10], [TT], [12]. This paper provides a first investigation 
of the rotopulsating orbits in §^ and H^. 

4. Results 

The remaining part of this paper is structured as follows. We first introduce 
in Section 5 the equations of motion in S'^ and M.^ and obtain their seven first 
integrals: one for the energy and six for the total angular momentum. In Section 
6, we define the rotopulsating orbits (or rotopulsators) and classify them into five 
types, which we call positive elliptic, positive elliptic-elliptic, negative elliptic, 
negative hyperbolic, and negative elliptic-hyperbolic, depending on the manifold 
in which they occur and on the nature and number of their rotations. 

The main results occur in Sections 7 through 13. Sections 7, 8, 10, 11, and 12 
have the same structure, each providing a criterion for the existence of the type of 
rotopulsator it analyzes, giving the corresponding conservation laws, and proving 
existence and uniqueness results for classes of such orbits in the case N = 3. Thus 
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Section 7 deals with positive elliptic, Section 8 with positive elliptic-elliptic, Sec- 
tion 10 with negative elliptic. Section 11 with negative hyperbolic, and Section 12 
with negative elliptic-hyperbolic rotopulsators. The examples provided in each of 
these sections are of Lagrangian type (i.e. the bodes lie at the vertices of rotating 
equilateral triangles) for the positive elliptic, positive elliptic-elliptic, and nega- 
tive elliptic rotopulsators, and of Eulerian type (i.e. the bodies lie on the same 
rotating geodesic) for the negative hyperbolic and negative elliptic-hyperbolic ro- 
topulsators. 

Sections 10 and 13 offer a theorem each, whose object is to describe the qual- 
itative behaviour of rotopulsators in and H^, respectively. The first theorem 
shows that, for rotopulsators of S^, for any foliation of the sphere with Clifford 
tori, none of the bodies can stay confined to some Clifford torus, but each body 
passes through a continuum of such surfaces. The second theorem proves a similar 
result for H^, where the Clifford tori are replaced by hyperbolic cylinders. 

An interesting finding is that of a class of Lagrangian orbits (see Subsection 8.3), 
which behave like relative equilibria, in the sense that they maintain constant 
mutual distances during the motion, but cannot be generated from any single 
element of the underlying subgroup 5*0(2) x 5*0(2) of the Lie group 50(4). We 
call these rotopulsating orbits nonstandard relative equilibria. 



Consider point particles (bodies) of masses ■mi,...,mN > moving in 
§^ (thought as embedded in the Euclidean space M'') or H'^ (embedded in the 
Minkowski space M^'-*^), where 



In previous work, we derived the equations of motion of the curved A^-body prob- 
lem using constrained Lagrangian dynamics and showed that by suitable coordi- 
nate and time-rescaling transformations the study of the problem can be reduced 
to §3 and [6], [7j. 

The configuration of the system is described by the 4A^-dimensional vector 



where = {wi, Xi, yi,Zi),i = 1, . . . , A^, denote the position vectors of the bodies. 
The equations of motion are given by the second-order system 



5. Equations of motion 



If 



= {{w, X, y, z) \ + x'^ + y'^ + 
{(w, X, y^z^ \ w'^ + x^ + y"^ — z^ = - 



= 1}, 

I, z> 0}. 



q = (qi,---,q7v), 



N 



[a-a(qi0q,)2]3/2 



(1) 



o"(qi qi)qi, « = i, . . . 



with constraints 



(2) 



qi qi = qi qj = 0, i = l,...,N. 
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Here © is the standard inner product ■ of signature (+, +, +, +) in or the 
Lorentz inner product □ of signature (+, +, +, — ) in H^, and 



a 



+1, in 
-1, in 



denotes the signum function. Equations ([T]) can be viewed as a 6A^-dimensional 
first-order system of ordinary differential equations. The force acting on each 
body has a gravitational component (the above sum) and a term (involving the 
velocities) that corresponds to the constraints. 

As a consequence of Noether's theorem, system ([T]) has the scalar integral of 
energy, 

T(q,q)-f/(q) = /i, 

where 

(Tmjmjqi © q^ 

^y'i)- 2^ f^_^(q 0g )213/2 
l<i<j<N L J 

is the force function {—U representing the potential), which stems from the cotan- 
gent of the distance. 



1 ^ 



i=l 

is the kinetic energy, with h representing an integration constant. System ([T]) also 
has the 6-dimensional integral of the total angular momentum, 

N 

mcii A qi = c, 

i=l 

where A is the wedge product and c = {c^x, c^y, c^z, Cxy, c^z, Cyz) denotes an inte- 
gration vector, each component measuring the rotation of the system about the 
origin of the frame relative to the plane corresponding to the bottom indices. 
For positive curvature, i.e. in we can explicitly write the equations of motion 

as 

(3) = E '^u'^'i ^"^^ '^i?^ - (q- • q-)q-' • = ^ = 1,2,. ..,iv, 

^ [1 - (q, ■ qj)2]^/2 



while for negative curvature, i.e. in H^, they take the form 



N 



4 qi = > Vn — V2 — Tiv^ + q^ ^ q* q^' q* h qi = -1, i = i,2,...,N. 

^ (qi □ qj)2 - If 



6 



Florin Diacu and Shima Kordlou 



Systems ([3]) and Q are the forms of the equations of motion we will mainly refer 
to in the remainder of this paper. 

6. Classification 

In this section we define several types of rotopulsating orbits of the curved A^- 
body problem and explain the choices we make to classify them. The rotopulsating 
orbit extends the Euclidean homographic solution to spaces of nonzero constant 
curvature, which actually reduce to the unit 3-sphere S^, in the positive case, 
and to the unit hyperbolic sphere, H^, for negative curvature. In two previous 
papers, we introduced this concept in the 2-dimensional case and kept using the 
name "homographic" for it, [5], [9]. Our idea was that the configurations we stud- 
ied (mostly polygons) remained homographic if viewed in the ambient Euclidean 
space. But it seems more natural to regard configurations in intrinsic terms, the 
more so when we move from two to three dimensions. 

In §^,§^,]HI^, and H^, however, the concept of similarity, which corresponds to 
the adjective homographic, makes little sense, since, for instance, the only similar 
triangles are the congruent ones. Consequently, to extend the concept of homo- 
graphic orbit to spaces of constant curvature, we need, on one hand, to capture 
somehow the dilation/contraction aspect as well as the rotation component of the 
solutions, and, on the other hand, to recover the original definition when the cur- 
vature tends to zero. We therefore introduce here a new adjective, rotopulsating, 
which suggests both these features of the orbit without implying similarity of the 
configuration. For rotopulsating orbits will be also use the noun rotopulsator. 

The definitions we provide below follow naturally from the concept of relative 
equilibrium of the curved A^-body problem, defined in [B], [7]. We introduced there 
various kinds of relative equilibria in terms of the isometric rotation groups of §^ 
and H^. The rotopulsating orbits we will further define differ from the relative 
equilibria through nonuniform rotations and, in general, through nonconstant 
mutual distances, as we will further see. 

Definition 1 (Positive elliptic rotopulsators). A solution of system ([s]) m 
is called a positive elliptic rotopulsator if it is of the form 

q = (qi,q2, • • • ,qiv), q^ = (wi,^^,?/^,^^), z = 1,2, . . . , A^, 

Wi = ri{t)cos[a{t) + ai], a;^ = ri(t) sin[a(t) + a^], yi = yi{t), Zi = Zi{t), 

where ai, i = 1,2, . . . , N , are constants, a is a nonconstant function, ri,yi, and 
Zi satisfy the conditions 

0<ri<l, -l<yi,Zi<l, and + y^ + z^ = 1, i = l,2,...,N, 

and Cyz = 0. If r is constant and a(t) = at, with a a nonzero constant, then the 
solution is called a positive elliptic relative equilibrium. 
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Remark 1. The conditions that a is nonconstant and Cyz = are imposed to 
ensure that the system has an elhptic rotation relative to the wa;-plane, but no 
eUiptic rotation relative to the yz-plane. Rotations relative to other base planes 
may occur. 

Definition 2 (Positive elliptic-elliptic rotopulsators). A solution of system 
([3]) in is called a positive elliptic- elliptic rotopulsator if it is of the form 

(6) q = (qi, q2, . . . , qAr), q» = (wj, Xj, y^, Zi), i = 1, 2, . . . , iV, 

Wi = Tiit) cos[a;(t) + Oj], Xi = Tiit) sin[a;(t) + Oj], 

Vi = Pi(t) cos[/3(t) + bi], Zi = pi(t) sm[f3(t) + bi], 

where ai,bi, i = 1,2, . . . , N , are constants, a and (3 are nonconstant functions, 
and Ti and pi satisfy the conditions 

0<ri,pi<l and r'f + p'^^ = 1, i = l,2,...,N. 

If r and p are constant and a{t) = at, = pt, with a, (3 nonzero constants, 
then the solution is called a positive elliptic- elliptic relative equilibrium. 

Remark 2. The condition that a and /3 are nonconstant is imposed to ensure 
that the system has two elliptic rotations, one relative to the wx-plane and the 
other relative to the yz-plane. Rotations relative to the other base planes may 
occur. 

Definition 3 (Negative elliptic rotopulsators). A solution of system ^ in 
is called a negative elliptic rotopulsator if it is of the form 

q = (qi,q2,...,q7v), qi = {wi,Xi,yi,Zi), i = l,2,...,N, 

Wi = ri{t)cos[a{t) ai], = ri(t) sin[a(t) + a^], yi = yi{t), Zi = Zi{t), 

where ai, i = 1,2, . . . , N , are constants, a is a nonconstant function, ri,yi, and 
Zi satisfy the conditions 

Zi>l and r^ + y^-z^ = -l, i = l,2,...,N, 

and Cyz = 0. If r is constant and a(t) = at, with a a nonzero constant, then the 
solution is called a negative elliptic relative equilibrium. 

Remark 3. The conditions that a is nonconstant and Cyz = are imposed to 
ensure that the system has an elliptic rotation relative to the wx-plane, but no 
hyperbolic rotation relative to the yz-plane. Rotations relative to other base 
planes may occur. 

Definition 4 (Negative hyperbolic rotopulsators). A solution of system Q 
in M.^ is called a negative hyperbolic rotopulsator if it is of the form 

q = (qi,q2, • • • ,qAf), qi = iwi,Xi,yi,Zi), i = l,2,...,N, 

Wi = Wiif), Xi = Xiit), yi = r]i{t) sinh[/3(t) + bi], Zi = r]i{t) cosh[/3(t) + bi]. 
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where bi, i = 1,2, . . . , N , are constants, (3 is a nonconstant function, Wi,Xi,Zi, 
and rji satisfy the conditions 

Zi>l and vol + x1 — rjl = —1, i = l,2,...,N, 

and Cwx = 0. If rj is constant and (3{t) = f3t, with (3 a nonzero constant, then the 
solution is called a negative hyperbolic relative equilibrium. 

Remark 4. The conditions that /3 is nonconstant and c^x = are imposed to 
ensure that the system has a hyperbohc rotation relative to the y^-plane, but no 
elhptic rotation relative to the wx-plane. Rotations relative to other base planes 
may occur. 

Definition 5 (Negative elliptic- hyperbolic rotopulsators). A solution of 
system Q in is called a negative elliptic-hyperbolic rotopulsator if it is of the 
form 

q = (qi,q2, • • • ,q7v), qi = {wi,Xi,yi,Zi), i = l,2,...,N, 

(9) Wi = Tiif) cos[a(t) + aj], Xi = Viit) sin[a(t) + aj], 

yi = r]i(t) sinh[/3(t) + bi], Zi = r]i(t) cosh[/3(t) + bi], 

where ai,bi, i = 1,2, . . . , N , are constants, a and /3 are nonconstant functions, 
whereas ri,rii, and zi satisfy the conditions 

Zi > 1 and - r]^ = -1, i = l,2,..., N. 

When r and rj are constant and a{t) = at, P{t) = (3t, with a, (3 nonzero constants, 
then the solution is called a negative elliptic-hyperbolic relative equilibrium. 

Remark 5. The conditions for a and /3 are imposed to ensure that the system 
has an elliptic rotation relative to the wx-plane and a hyperbolic rotation relative 
to the ?/2;-plane. Rotations relative to other base planes may occur. 

Remark 6. There are alternative ways to define the above classes of rotopulsating 
solutions of the curved A^-body problem (including those of considering positive 
elliptic orbits as particular positive elliptic-elliptic orbits, and negative elliptic 
and negative hyperbolic orbits as particular negative elhptic-hyperbolic orbits), 
but our choice is more convenient for computations and will help us emphasize, 
unambiguously, certain properties of all these solutions. 

Before starting to study these orbits, let us remark that, in spaces of nonzero 
constant curvature, the main difference between relative equilibria and rotopul- 
sators is the same we find in Euclidean space between relative equilibria and 
homographic orbits that are not relative equilibria. Namely, in the former class of 
orbits the mutual distances are constant in time, whereas in the latter class, the 
geometric configuration of the particle system changes in size, in general. This is 
the reason why we used in [6] and [7] the isometric rotation groups of and 
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to define the relative equilibria, a move that led us naturally to five classes of such 
orbits. Allowing then the size of these configurations to vary, and implicitly mak- 
ing the rotation of the system nonuniform, we were led to the above definitions 
for rotopulsators. 

In 3-dimensional spaces of constant curvature, however, a new possibility oc- 
curs, which does not take place in ]R^,M^,S^, and IP. Namely, there could exist 
rotopulsating orbits, as defined above, which maintain constant mutual distances. 
In other words, such motions (if they exist, and we will show that they do) are 
relative equilibria, though they don't belong to the class of orbits obtained from 
isometric rotations, i.e. do not obey the definitions of relative equilibria given 
above, or in [6] and [7]. Therefore the following new definition becomes necessary. 

Definition 6. Any of the rotopulsators defined in (|6|, or Q are 

called nonstandard relative equilibria if they are not relative equilibria in the sense 
of Definitions [T| |2], |3], |4], or [5| although the mutual distances between the bodies 
remain constant during the motion. 

As we will see later, these orbits occur due to the fact that they cannot be 
generated from a single element of the Lie group 5*0(2) x S0{2), which occurs 
naturally when describing such solutions. Nevertheless, they may occur from an 
element generated by 5*0(4), as the more general theory of relative equilibria 
developed in the framework of geometric mechanics suggests, [T9] . 

7. Positive elliptic rotopulsators 

In this section we analyze the solutions given in Definition [T} We first introduce 
a criterion for finding them, then obtain the conservation laws, and finally prove 
the existence of a particular class of orbits, namely the positive elliptic Lagrangian 
rotopulsators of the 3-body problem in 

7.1. Criterion for positive elliptic rotopulsators. The following result pro- 
vides necessary and sufficient conditions for the existence of positive elliptic ro- 
topulsators in S^. 

Criterion 1. A solution candidate of the form ^ is a positive elliptic rotopulsator 
for system ^ if and only if Cyz = 0, 

(10) a 



where c is a nonzero constant, and the variables yi, Zi, i = 1,2, . . . , N, satisfy the 
system of 2N second-order differential equations 



(11) 



_ sr^N mj{zj-qijZ,) _ p 
i Z^j=l ("l— o2 )3/2 ^iZi, 



10 
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yf + - iviZi - Ziijif 



+ 



h ^3 



2 ' 



^ = 1,2, 



and, for any i,j G {1, 2, ... , A^}, qij is given by 

lij := Qi ■ qj = (1 - y1 - zl)^{^ - Z/j - z]Y^ cos(ai - a^) + yiy^ + ZiZj. 

Proof. Consider a solution candidate of the form ([s]). A straightforward compu- 
tation shows that, for any i, j G {1, 2, . . . , A^}, is of the form given in the above 
statement. Moreover, for any i = 1, 2, . . . , A^, we find that 

. . _ yf + zf - {yjZj - Zjijif + (1 - yf - zffa^ 
For all z = 1, 2, . . . , A^, each rj can be expressed in terms of yi and Zi to obtain 



;i - yl - z, 



2)5 
I I 



(ViZi - Ziijif - ijf - zf - (1 - yf - zj){yiyi + z^iJ^) 



:^-y 



2 -^^Nf 



Substituting a candidate solution of the form (|5| into system ([s]) and employing 
the above formulas, we obtain for the equations corresponding to yt and Zi that 



N 



(13) m 



T^jiyj - Qijyi) Wi + Zi - (yiZi - ZiyiY\lH 2 _ 2^ -2 

\^ yi z^ jyicx , 



-4) 



TV 



(14) 



'^ji.Zj - QijZi) [yf + Zi - {yiZi - z^y^y]zi 



^ fl- 



i-yf-zf 



- (1 - - zDziO^. 



For the equations corresponding to ivi and Xj, after some long computations that 
also use (13) and (14), we are led either to identities or to the equations 

(15) 



N 



ridi + 2rid 



m,ri sm(aj — Oj 



^ (1 



l,2,...,Ar. 



For every i = 1,2, . . . , N , we further multiply the ith. equation in (15) by miVi, 
add the resulting A^ equations, and notice that 

y~v y~v mirrijrirj sm(aj — ajj ^ 



:i-?^ 
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N 



N 



a 



miVi + 2d ^ mirif-i = 0, 



1=1 



1=1 



which has the solution 

a = 



where c is an integration constant. Consequently, equations ([13j) and (14) become 
(16) 



N 



(17) 



^jiVj - QijVi) [yf + - iViZi - Ziyif]yi c^{l - yf - zf)yi 



N 



2ij, 



^ji^j - Qij^i) [yf + - hji^i - Z^yif]z^ 



1 - Z/i - z, 



[Ej=imj{l-y]-z] 



2 ' 



2 ' 



i = 1,2, . . . , N. Since equations ( |16| ) and (17) lead to system (11), and the 
condition Cyz = occurs from Definition [T| the criterion is proved. □ 

7.2. Conservation laws for positive elliptic rotopulsating orbits. In addi- 
tion to Criterion [T| we would also like to obtain the conservation laws specific to 
positive elliptic rotopulsators. They follow by straightforward computations using 
the above proof, the integral of energy, and the six integrals of the total angular 
momentum. 

Proposition 1. If system (|3| has a solution of the form then the following 
expressions are constant: 
— energy, 

^ rriilyf + zf - {yiZi - Ziyi)"^] 



i=l 



2(1 -yf-zj 



+ 



2ELmil-yf- 



^ ^ AT I -L 



mirrijqij 



.2\i 



l<i<j<N 



— total angular momentum relative to the wx-plane, 
(19) c. 



where c is the constant in the expression (15) of a; 



12 Florin Diacu and Shima Kordlou 

— total angular momentum relative to the wy -plane: 

rrii a-Vi - Zi)^yi+ -j 



(20) 



N 

E 

1=1 



cos(q; + ai) 



N 



^mi(l - - Zi)^yiSm{a + ai); 



2 

i ! i=l 



total angular momentum relative to the wz -plane, 



(21) 



N 

E 

i=l 



rrii 



zf)-^Zi + 



iViVi + ZiZi)zi 



cos(a + tti) 



+ 



E7=i"^j(1-1/j 



N 

^mi{l-y1 - Zi)^Zi sm{a + a^); 



^3 ] ' j=l 

total angular momentum relative to the xy-plane, 

N 



--xy 



(22) 



E 

1=1 



rrii 



sin(Q; + tti) 



N 



^mi(l - yl - zlY^yiCO&{a + a^); 



J2j=imj{l-y]-z]) 
total angular momentum relative to the xz -plane: 

{yiyi + ZiZi)zi 



(23) 



AT 

E 

1=1 



sva{a + Cj) 



AT 



(24) 



E,=i^i(i-y|-^|) i=i 

total angular momentum relative to the yz-plane: 

c =0 



The above result could be also used to prove the nonexistence of some candi- 
dates for positive eUiptic rotopulsators, by showing that at least one of the above 
conservation laws is violated. 

7.3. Positive elliptic Lagrangian rotopulsators. We further provide a class 
of specific examples of positive elliptic rotopulsators of the curved 3-body prob- 
lem, namely Lagrangian orbits, i.e. bodies that lie at the vertices of a rotating 
equilateral triangle in that changes size, which means that it fails to be similar 
to itself but has congruent sides at every instant in time. As we will see, these 
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systems rotate relative to the wx-plane, but have no rotations with respect to the 
other base planes. 

Consider three equal masses, mi = m2 = ='■ m, and a candidate solution of 
the form 

(25) q = (qi,q2,q3), qi = iwi,Xi,yi,Zi), i = 1,2,3, 

Wi = r(t) cosa(t), Xi(t) = r(t) sina;(t), yi = y(t), Zi(t) = z(t), 
W2 = r(t) cos[a(t) + 27r/3], X2{t) = r{t) sm[a{t) + 27r/3],i/2 = y{t), Z2{t) = z{t), 
W3 = r(t)cos[a(t) +47r/3], x^it) = r(t) sin[a(t) + 47r/3], ?/3 = y{t), Zs{t) = z{t). 

With the help of Criterion JTJ we will further show that, in general, these are 

indeed solutions of system (|3j. 

Proposition 2. Consider the curved 3-body problem in given by system ([s]). 
Then, except for a negligible set of solutions formed by positive elliptic Lagrangian 
relative equilibria, every candidate solution of the form (25) is a positive elliptic 
Lagrangian rotopulsator, which rotates only relative to the wx-plane. 



Proof. Let us consider a candidate solution of the form (25 ). Then, using Criterion 
[T| straightforward computations show that 



Qu — 121 — lis — 131 — Q23 — Q32 



3y^ + 3z' 



a 



2 ' '~ 3m(l 
and that the variables y and z must satisfy the equations 

y = F{y,z,y,z)y 
z = F{y,z,y,z)z, 



y2 _ ~;2j 



(26) 
where 



F{y,z,y,z) 



8m 



V3{1 



z^)-2{l + 3y^ + 3z^) 
{yz - yzf 



9m2(l — y"^ 



y"^ + z^ 



From (26), we can conclude that yz = yz, which implies that 

yz — zy = k (constant). 



But, since from (24) we have that 3m(yz — zy) = Cyz, it follows that k = Cyz/3m. 
As Cyz = 0, we must take k = 0, so yz — zy = 0, and therefore J^^ = if 2; does 
not take zero values, so y{t) = 'yz{t), where 7 is a constant. Moreover, since 

sina + sin(a + 27r/3) + sin(a + 47r/3) = cos a + cos(q; + 27r/3) + cos(a + 47r/3) = 0, 

V? = 0, so the 



it follows from (|20j), ([21]), (l22j), and (l23j) that c 
triangle has no rotation relative to the planes wy, wz, xy, xz, and yz. 



wy 
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Notice that the energy relation (18) takes the form 



3m[|/2 + i2 _ - ijzf] \/3m2(3|/2 + Sz^ - 1) 

2(1 - y2 _ ^2) +6r„(l - y2 _ ^2) - _ ^2 _ ^2)1(1 + 3^2 + 3^2) | 

which imphes that F can be written as 

2m[5 - 9(y2 + ^2)2] 2/1 



i^(2/,^) 



v^(l-|/2_;22)|(i + 3^2 + 3^2)| 3m 



If we denote 5 = 7+1 > 1, system (26) reduces to the family of first order 
systems 



(27) 



2m(5-9(52z4) 



3m 



.V^(l-<522) 2 (1+3(5^2)2 

Standard results of the theory of ordinary differential equations can now be 



applied to system (27) to prove the existence and uniqueness of analytic positive 
elliptic Lagrangian rotopulsators, for admissible initial conditions. To show that 
only a negligible set of such orbits is given by relative equilibria, notice that finding 



the fixed points in the family of equations given by system (27), which depends 
on parameters, reduces to obtaining the zeroes of a polynomial equation. This 
remark completes the proof. □ 

7.4. The nature of the solutions. Although we proved the existence of positive 
elliptic rotopulsators, it would be still interesting to learn more about their nature. 



For this, we will first find the fixed points for the vector field of system (27) 
different from the obvious one, (2;, v) = (0, 0), which is common to all equations in 
the family. This problem reduces to finding the zeroes of the family of polynomials 

P{z) = 27(9m^ + /i2)5V - 18(15m^ + h^)5^z^ - 8h^5z^ + 75m^ - h^. 

According to Descartes's rule of signs, we have to distinguish between two cases: 

(i) \h\ < 5\/3m'^, when, for every fixed values of the parameters, P has either 
two positive or no positive root at all; 

(ii) 1^1 > by/Sm^, when, for every fixed values of the parameters, P has exactly 
one positive root. 

Case (ii) always leads to one fixed point since the unique positive root, zq := 
zo{m, h, e), is always such that \zo\ < 1. To show this fact, let ( = 6z'^. Then we 
can associate to P the family of polynomials 

P(C) = 27(W + h'^X^ - 18(15m^ + h'^)C - 8h\ + 75m^ - /il 

Since 5 > 1, to prove that |2;o| < 1 it is enough show that |Co| := If^^ol < 1- For 
this, notice that P(l) = A8m^ > 0. As P(0) < in case (ii), it follows that 
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< Co < 1- Then the corresponding eigenvalues Ai,2 are given by the equation 



2h 

X+ Wizo,m,e) = 0, 

Sm 

where W{zo,m,e) is a finite number for every fixed value of the parameters. In- 
dependently of the values of W{zQ,m,e), the eigenvalues show that z varies for 
every orbit that is not a fixed point. Similar conclusions can be drawn in case (i). 



Numerical experiments suggest that all the other orbits of system (27) are peri- 



odic, except for two homoclinic orbits. Since the periods of a and z don't usually 
match, the periodic orbits generate quasiperiodic positive elliptic Lagrangian or- 
bits, except for a negligible class, given by periodic positive elliptic Lagrangian 
orbits (see Fig. [T]) . 




Figure 1. The fiow of system (27) for m = l,Cyz = 0,6 = l,h = 5. 



Remark 7. When y or z is a nonzero constant, the motion takes place on a 
2-dimensional nongreat sphere. When y = or z = 0, the motion takes place on 
a 2-dimensional great sphere, i.e. on a type of Lagrangian orbit for which we 
gave a complete classification in [9]. 

8. Positive elliptic-elliptic rotopulsators 

In this section we analyze the solutions given in Definition |2} We first prove 
a criterion for finding such solutions, then obtain the conservation laws, and fi- 
nally discuss a particular class of examples, namely the positive elliptic-elliptic 
Lagrangian rotopulsators of the 3-body problem in S^. 

8.1. Criterion for positive elliptic-elliptic rotopulsators. The following re- 
sult provides necessary and sufficient conditions for the existence of positive 
elliptic-elliptic rotopulsators in S^. 
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Criterion 2. A solution candidate of the form (|6j) zs a positive elliptic- elliptic 
rotopulsator for system (|3|) if and only if 



(28) 



a 



Cl 



EN 2 
i=l ^i^f 



/3 



C2 



with ci,C2,M = YliLi ™i nonzero constants, and the variables ri, r2, . . . , satisfy 
the N second-order differential equations 



Ti = ri{l - r1) 



(29) 



- (Ei=i "nfiirff (M - Y.i=i rmrl 



■ 2 

rjri 
1-r? 



N 



mj[rj{l — rf) cos(aj — aj) — rj(l — rf)i{l — r^)2 cos(6j — bj)] 



(1-4) 



where, for any i,j G {1, 2, ... , A^} wzi/i i ^ j , we denoted 



(30) 



qi ■ qj- = r^rj cos(ai - aj) + (1 - r,^) ^ (1 - rj) 2 cos(6i - bj). 



Proof. Consider a candidate solution of the form (|6| for system ^ . By expressing 
eacli Pi in terms of rj, z = 1, 2, . . . , N, we obtain that 



Pi = (l-^r)^ P 



'1 — r, ) 2 



Pi 



r,^ + ri(l -rf)fj 
(l-r?)-2 ' 



eij takes the form (30), and 



q. ■ q. = + r^a' + + (1 - rf)/?^ 

1 — r,- 



Substituting a solution candidate of the form ^ into system ([3]), and using the 
above formulas, we obtain for the equations corresponding to Wt and Xj the equa- 
tions 



l-ri 



(31) 



N 



mj[rj{l — r'j) cos(ai — a^) — rj(l — rj^)2(l — r^)2 cos(6j — bj)] 



(1-4) 



(32) r,a + 2f,d = -V ' ] \ '\ z = 1, 2, . . . , iV, 

i:i (1-4)' 
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respectively, whereas for the equations corresponding to jji/zi, we find equations 
(31) again as well as the equations 

N 

(33) + 2r,/3 



Em,(l - r^)2 sin(6i - bj 



^ = 1,2, 



AT. 



We can solve equations (32) the same way we solved equations (15) and obtain 

Cl 



a 



i 



where Ci is an integration constant. To solve equations (33), we proceed similarly, 
with the only change that for each i = 1,2, . . . , N, the corresponding equation 
gets multiplied by mj(l — r^)? instead of m^rj, to obtain after addition that 

C2 



/3 



where M = YliLi^i ^2 is an integration constant. Then equations (31) take 
the form ( 29 ) , a remark that completes the proof. □ 



Remark 8. From (28), we can conclude that a and (3 are not independent of 
each other, but connected by the relationship 

(34) 



a (3 



written under the assumption that a and P are not constant. In particular, if a 
and P differ only by an additive constant, they are linear functions of time, i.e. 



(35) 



a = /3 



Cl + C2 

M 



Then system (29) becomes 
(36) 



. 2 



N 



■mj[rj{l — rf) cos(aj — aj) — rj(l — rf) 2 (1 — rj) 2 cos(6i — hj)] 



;i-4) 



1,2 N. 



8.2. Conservation laws for positive elliptic-elliptic rotopulsators. In ad- 
dition to Criterion [2| we would also like to obtain the conservation laws specific to 
positive elliptic-elliptic rotopulsators. These laws follow by straightforward com- 
putations using the above proof, the integral of energy, and the six integrals of 
the total angular momentum. 
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Proposition 3. If system (|3| has a solution of the form then the following 
expressions are constant: 
— energy: 



N 



(37) h = 



E 



rriimjeij 



— total angular momentum relative to the wx -plane, 
(38) = Ci, 



where ci is the constant in the expression (28) of a; 
— total angular momentum relative to the wy-plane, 

N 



1=1 



rj(l — r1) 2 (d + /3) sin(Q; — (5 + — hi) 



1 ' I 

(39) +rj(l — rj^)2 (d — (3) sm{a + [3 + + hi) ^ — ^ cos(a; — (3 + Qi — hi) 

(1 -r2)2 

— cos(a + l3 + Ui + hi) 



(1 - 

total angular momentum relative to the wz-plane, 



1 ^ 



nii 



i=l 



rj(l — r^) 2 (d + /3) cos(a + f3 + ai + hi 



(40) +ri{l - rf)^a - /3) cos(a - f3 + ai - hi) + - — sin(a - f3 + ai - k) 

(l-rf)2 



'-ir^ sin(a + (3 + ai + hi) 

(l-rf)2 

total angular momentum relative to the xy-plane, 

N 



xy 



4e 



i=l 



rrii 



Tiil - r1) 2 (d + /3) cos(a; - /3 + - 6^ 



(41) +ri(l - r,^)^(d - /3) cos(a + (3 + ai + hi) + '—^ sin(a - /3 + a^, - 6. 

(l-r/)2 



Y sin(a + /3 + ttj + 6j) 
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total angular momentum relative to the xz-plane, 



1 ^ 

E 



rrii 



i=l 



rj(l — rf) 2 (d + /3) sin(a — (3 + — bi 



(42) -r,(l - rf)^a - /3) sin(a + /3 + + 6,) 



1 — rf)2 



Y cos(a — (3 + tti — bi) 



"—T cos{a + (3 + ai + bi 

total angular momentum relative to the yz-plane, 



(43) 



C2, 



where C2 ^ is the constant in the expression (28) of j3 



The above result could be also used to prove the nonexistence of some candi- 
dates for positive elliptic rotopulsating orbits, by showing that at least one of the 
above conservation laws is violated. 

8.3. Positive elliptic-elliptic Lagrangian rotopulsators. We further provide 
a class of specific examples of positive elliptic-elliptic rotopulsators of the curved 
3-body problem, namely Lagrangian orbits, i.e. bodies that lie at the vertices of a 
rotating equilateral triangle in S'^. These systems rotate with respect to the wx- 
and the y^-plane, but have no rotations relative to the other base planes. These 
solutions provide the first examples of nonstandard relative equilibria, as stated 
in Definition 6, i.e. every such equilateral triangle remains congruent to itself, in 
spite of the fact that its orbit cannot be generated by the action of any single 
element of the underlying Lie group 5*0(2) x S0{2), which occurs naturally from 
the way we fixed the coordinate axes, although it may be generated by an element 
of 50(4). 

Consider three equal masses, mi = m2 = rn^ =: m, and a candidate solution of 
the form 

(44) q = (qi,q2,q3), qi = iwi,Xi,yi,Zi), 2 = 1,2,3, 

wi = r{t) COS a(t), xi(t) = r(t) sin a(t), yi = p(t) cos (3 (t), zi(t) = p(t) sin (3 (t), 
W2 = r{t) cos[a{t) + 2n/3], Xiit) = r{t) sin[a{t) + 27r/3], 
y2 = p{t) cos[(3{t) + 27r/3], Z2{t) = p{t) sin[(3{t) + 27r/3], 
W3 = r(t) cos[a(t) + 47r/3], X3(t) = r(t) sin[a{t) + 47r/3], 
7/3 = p{t) cos[/3(t) + 477/3], z,Xt) = p{t) sin[/3(t) + 47r/3], 

with a and [3 nonconstant functions and + p'^ = 1. We can now prove the 
following result. 
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Proposition 4. Consider the curved 3-body problem in §^ given by system ([s]). 
Then, except for a negligible set of solutions formed by positive elliptic- elliptic 



Lagrangian relative equilibria, every candidate solution of the form (44) is a posi- 
tive elliptic- elliptic Lagrangian rotopulsator and, at the same time, a nonstandard 
relative equilibrium. 

Proof. To prove this result, let us first apply Criterion |2] to a candidate solution 
of the form (44) for system ([3]). Straightforward computations show that 



ei2 — ^21 — ei3 — ^31 — ^23 — ^32 — ^1/2) 

which means that the sides of the equilateral triangle don't vary in time, so if this 
solution candidate proves to be a positive elliptic-elhptic Lagrangian rotopulsator, 
then it is also a nonstandard relative equilibrium. 
Notice further that 

(45) « = ;7^, P 



Smr"^ 3m(l — r^) 



and that system (29) reduces to the family of first-order systems 
(46) 



^ 9m^r'' 9m^(l— r^) 



As in the proof of Proposition [2} the existence and uniqueness of analytic pos- 
itive elliptic-elliptic rotopulsators for admissible initial data follows. □ 

8.4. The nature of the solutions. To get some insight into the nature of these 



solutions, let us first find the fixed points of system (46), which occur for 

t 2 2\4 o22|2 n 

(Ci — Cajr — 2cir -|- q = 0. 



For Cl 7^ C2, there are at most two fixed points, ri = ^ ^^^^^ and r2 — y 

only one of which is between and 1, whereas for ci = C2 there is a single fixed 
point, namely tq = 1/a/2. So in both cases the fixed point is unique. 



Numerical experiments suggest that the phase space picture of system (46 ) looks 
like in Fig. ([2]), with r periodic, a fact in agreement with the Lie theory applied 
to the group 5*0(4). Consequently, the positive elliptic-elliptic Lagrangian orbits 
are expected to be not periodic, since the periods of a and (3 differ, in general, 
except for a negligible set corresponding to periodic orbits. 



Remark 9. From relations (45) we can see how the angular velocities a and f3 
vary relative to r: when r is close to 0, |d| is large, while is small, and the 
other way around when r is close to 1. Moreover, the expressions of a and P are 



consistent with relation (34) 
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Figure 2. The flow of system (46) for m = 1/3 and ci = C2 = 1 



For Cl 7^ C2, tlie flow looks qualitatively similar. 

Remark 10. If we initially assume the masses to be distinct, but the bodies to 
correspond to the same function r, it follows that the solutions don't exist, which 
means that the masses must be always equal. If we additionally take distinct 
functions ri,r2,r3, then 

r.Tj + (1 - r2)V2(l _ ^2)1/2 _ ^.)2 ^ _ ^2)1/2 _ _ 2)l/2]2 _ 3 



''2 4 ' 

which implies that the triangle is not necessarily equilateral, so the orbit may not 
be Lagrangian at all. 



Remark 11. The energy relation (37) takes the form 

3mr^ 



h 



1 



2(1 — r^) Qm? 
so the energy constant, /i, is always positive. 

Remark 12. When d = /3, we have = 



constant, if it exists, so 



a 



/3 



C1+C2 ' 

Cl + C2 

3m 



which means that r must be 



These facts imply that the orbit becomes a standard relative equilibrium, not a 



rotopulsator, and the integrals (37), (39), (40), (41), and (42) are constant, as 
expected. 

It is natural to ask whether there exist other nonstandard relative equilib- 
ria beyond the positive elliptic-elliptic Lagrangian rotopulsators described in this 



section. A close look at expression (30) leads to the following conclusion. 



Remark 13. If for for any positive elliptic-elliptic rotopulsator, i.e. any solution 
of the form ([g]) of system (|3]), we have = for all j G {1, 2, ... , N} and = 
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bi,i = 1, 2, . . . , A^, then this solution is also a nonstandard relative equilibrium. 
In particular, positive elliptic-elliptic rotopulsators given by regular polygons can 
satisfy this property. 

9. Qualitative behaviour of the positive rotopulsators 

In this section we will prove a result that describes the qualitative behaviour of 
rotopulsators in For this purpose, we first briefly introduce an object familiar 
to geometric topologists and present some of its properties. 

9.1. Clifford tori. The 2-dimensional manifold defined by 

(47) Trp:= {{w,x,y,z) eR'^ \ + = 1, 0<9,(f)<27i}, 

where w = rcos6,x = rsin6,y = pcos0, and z = psin0, with r, p > 0, is called 
a Clifford torus, and it has zero Gaussian curvature. Since the distance from to 
every point of T^p is 1, it follows that Clifford tori are included in S^. When r 
(and, consequently, p) takes all the values between and 1, the family of Clifford 
tori such defined foliates Each Clifford torus splits §^ into two solid tori and 
forms the boundary between them. The two solid tori are congruent only when 
r = p= l/^/2. 

We have previously shown that positive relative equilibria rotate on Clifford 
tori, in], [7j. We will next prove that, at every moment in time, a rotopulsator 
passes through a different Clifford torus of any given foliation of In other 
words, rotopulsators cannot be generated by an element of the underlying Lie 
group ^0(2) X S0{2). 

9.2. Positive rotopulsators. We can now state and prove the following result. 

Theorem 1. Consider a positive elliptic or positive elliptic- elliptic rotopulsator 
of the curved N-body in , for which a coordinate system of the ambient space 
is fixed. Then, for the corresponding foliation of given by Clifford tori, 
the trajectory of each body cannot be contained in any single Clifford torus, but 
intersects a continuum of Clifford tori. For positive elliptic and positive elliptic- 
elliptic relative equilibria, however, each body moves on some Clifford torus of the 
foliation for all time. 

Proof. The qualitative behaviour of positive relative equilibria in S^, including 
properties specific to each type of orbit, was previously described and proven to 
take place, [6], [7j, so we will further focus on positive rotopulsators. Since a 
reference frame was a priori fixed in M"^, the natural foliation of is given by the 
family (Trp)o<r,p,<i of Clifford tori, with T^p defined in (47). 



But the positive elliptic-elliptic rotopulsators are given by equations (|2]), with 
r and p nonconstant functions satisfying the relation + p^ = 1. Each Clifford 
torus in the foliation, however, corresponds to a constant r and a constant p that 
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also satisfy the relation + p"^ = 1. Consequently, at every moment in time, 
each particle passes through a different Clifford torus, thus travelling through a 
continuum of such objects during its motion. 

In the case of positive elliptic rotopulsators, the rotation takes place only with 
respect to the variables w and x, so we can treat this situation as a particular form 
of positive elliptic-elliptic rotopulsators, for which the result was proved above. 
This remark completes the proof. □ 



Remark 14. The nonstandard relative equilibria (44) behave like the positive 
elliptic-elliptic Lagrangian rotopulsators, not like the positive elliptic-elliptic La- 
grangian relative equilibria, in the sense that they pass through a continuum of 
Clifford tori instead of moving all along a single Clifford torus, as standard rela- 
tive equilibria do. But like the relative equilibria, they maintain constant mutual 
distances during the motion. 

10. Negative elliptic rotopulsators 

In this section we analyze the solutions given in Definition [3j We first prove a 
criterion for finding such solutions, then obtain the conservation laws, and finally 
discuss a particular class of examples, namely the negative elliptic Lagrangian 
rotopulsators of the 3-body problem in H^. 

10.1. Criterion for negative elliptic rotopulsators. The following result pro- 
vides necessary and sufficient conditions for the existence of negative elliptic ro- 
topulsators in H^. 

Criterion 3. A solution candidate of the form ^ is a positive elliptic rotopulsator 
for system ^ if and only if Cyz = 0, 

(48) a = —j^ , 

E,=i"^j(^|-2/|-l) 

where b is a nonzero constant, and the variables yi, Zi, i = 1,2, . . . , N, satisfy the 
system of 2N second-order differential equations 



(49) 



where 



(50) 



s-^N mj{zj+iJ.ijZi) I p 

(n'2 —1)3/2 ' -'^i'^i' 



[{yiZ, - Ziijif + zt - ijt] , b^{zt-yf 



^f-yf-^ [EUrn,{z]-y]-l)f 
i = 1,2, . . . , N, and, for any i,j G {1, 2, ... , A^}, /ijj is given by 

l^ij := qi ■ qj = {Zi - Vi - ^)Hz] - - 1)^ cos(ai - aj) + yiyj - ZiZj. 
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Proof. Consider a solution candidate of the form ([T]) subject to the above initial 
conditions. Then, for any i,j G {1, 2, ... , N}, we have 

fiij := Qi □ qj = {zf - y1 - l)^{z] - - 1)5 cos(a, - a^) + yiy^ - ZiZj, 
and, for any i = 1, 2, . . . , A^, we find that 

{ViZi - Zii/if + zf - ijf + {zf - yf - Ifa^ 



□ qi 



For all z = 1, 2, . . . , A^, each can be expressed in terms of yi and Zi to obtain 

ri = [Zi -yi -1)'^, ri = — -T, 

{zf-yt-iy^ 
.. ^ {zf - yf - l){z{Zi - yijji) + yf - zf - {yiZj - Zijjif 

{zf-yf-l)^ 

Substituting a solution of the form ([T]) into system Q and employing the above 
formulas, we obtain for the equations corresponding to iji and Zi that 



N 



N 



zf-yf-i 



(52) z, = Y^ m,(z, + iii^Zi) ^ [{y,z, z^if + 4^ yl\z^ ^ _ ^2 _ ^^^^^2 



whereas for the equations corresponding to Wi and Xj, after some long computa- 
tions that also use (51 ) and (52), we are led either to identities or to the equations 

(53) 



N 



EmiTi sm a, — a,- 
iul-l 



For every i = 1,2, . . . , N , we further multiply the ith equation in (53) by m^rj, 
add the resulting N equations, and notice that 

N N . , X 



mimiriTi sm[aj - ai) _ 

3 ~ u. 



(n^. - 1 

Thus we obtain the equation 

A 

a rriirf + 2a miTiTi = 



N 



1=1 



1=1 
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which has the solution 



a 



where b is an integration constant. Consequently, equations (51) and (52) become 
(54) 



N 



Vi 



mjiyj+HijVi) [{viZi - Ziijif + zf - yf]yi 9{zf -yf -l)y, 



(4-1) 



(55) 



N 



i:: (4-1)' 



+ 



mj(zj + HijZi) \{yiZi - Ziijif + zf - yf]zi b^{zf - yf - l)z 



zf-yf-i 



+ 



[Ef=i^jiz] - y] - 1)]^ 



i = 1,2, . . . , N. Since equations (54) and (55) are equivalent to system (49), and 
we have Cy^ = from the definition of negative elliptic rotopulsaating orbits, the 
criterion is proved. □ 

Remark 15. Notice that the constant c^x of the total angular momentum is 

N N 

Cwx = mi{wiXi - WiXi) = a mi{zf - y1 - 1) = 6, 

i=l i=l 

SO b describes the rotation of the particles relative to the wx-plane. 

10.2. Conservation laws for negative elliptic rotopulsating orbits. In ad- 
dition to Criterion [3} we would also like to obtain the conservation laws specific 
to negative elliptic rotopulsating orbits. They follow by straightforward compu- 
tations using the above proof, the integral of energy, and the six integrals of the 
total angular momentum. 

Proposition 5. If system ^ has a solution of the form then the following 
expressions are constant: 
— energy, 



(56) 



N 

i=l 
6^ 



rriiliyiZ, - Ziijif + zf - ijf] 



2(^^ -yf-i) 
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— total angular momentum relative to the wy -plane, 

iViVi - ZiZi)yi 



--wy 



(57) 



N 

1=1 



cosia + tti 



N 



total angular momentum relative to the wz -plane, 



(58) 



TV 
i=l 



{zt-yi-l)-^Zi + 



I . , {Viili - ZiZi)zi 



cos(q; + tti) 



N 



EN / 2 2 

i=imi[zi-yt 



— ^ mi{zl - y1 -if^Zi sin(a + a^); 

i=i 



total angular momentum relative to the xy-plane, 



N 



(59) 



{yjiji - ZiZ,)yi 



(z? 



sm{a + Oj) 



N 



^ mi{Zi - yl - l)^2/i cos(a + a^); 



Y.i=im{zl - yf - I) ^ 
total angular momentum relative to the xz -plane. 



(60) 



AT 



ITLi 



{z1-y1-l)-^Zi + 



{yjiji - ZiZi)zi 

iA-yl-^ 



sm(Q; + Oi) 



N 



m^(2;,^ - y,^ - 1)^2;, cos(q; + a^); 



Ei=i - yl - 1) 
— total angular momentum relative to the yz-plane, 

N 

(61) Cy^ = ^ mi{yiZi - yiZi). 

i=l 

10.3. Negative elliptic Lagrangian rotopulsators. We further provide a class 
of specific examples of negative elliptic rotopulsators of the curved 3-body prob- 
lem, namely Lagrangian orbits in H^. These systems rotate relative to the wx- 
plane, but have no rotations relative to other base planes. 

Consider three equal masses, mi = m2 = =: m, and a candidate solution of 
the form 



(62) 



Iqi,q2,q3j 



{wi, Xi, yi, Zi), i — 1,2, 3, 
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wi = r{t) cosa(t), xi(t) = r(t) sina(t), yi = y(t), zi(t) = z(t), 
W2 = r{t) cos[a(t) + 27r/3], X2{t) = r{t) sin[a(t) + 27r/3], ?/2 = y{t), Z2{t) = z{t), 
ws = r(t)cos[a(t) +477/3], Xsit) = r(t) sin[a(t) + 47r/3], 2/3 = y{t), z^{t) = z{t). 

With the help of Criterion |3} we can now show that, in general, these are indeed 
solutions of system Q. 

Proposition 6. Consider the curved 3-body problem in given by system Q. 
Then, except for a negligible set of solutions formed by negative elliptic Lagrangian 



relative equilibria, every candidate solution of the form (62) is a negative elliptic 
Lagrangian rotopulsator. 



Proof. Let us consider a candidate solution of the form (62 ). Then, using Criterion 
[3} straightforward computations show that 

3?/2 _ 3^2 ^ 1 ^ h 



/^12 — A'lS — /^23 



a 



2 ' 3m{z'^ 
and that the variables y and z satisfy the equations 

y = G{y,z,y,z)y 
J = G{y,z,y,z)z, 

where 



(63) 



G{y,z,y,z) 



{yz - yzf + z' 



y'-i 

8m 



+ 



Qm?(z'^ 



From (63), we can conclude that yz = yz, which implies that 

yz — zy = k (constant). 



But, since from (61) we have that 3m{yz — zy) = Cyz, it follows that k = Cyz/3m. 



Notice that the energy relation (\56v takes the form 



3m[{yz — yzY + z^ — y"^ 



+ 



^m\3y^ - 3z^ + 1) 



2(^2 -y^-l) 



6m{z'^ — y"^ 



1/2 -1)1(3^2 -3^2 + 1)1 



which implies that G can be written as 

^ 2/1 2m[5-9(y2 

^^y^'^ = ^ 

Since 



.2n2i 



y3(^ 



y' 



1)5(3^2 _ 3^2 



i)i 



sina + sin(a + 27r/3) + sin(a + 47r/3) = cos a + cos(q; + 27r/3) + cos(Q; + 47r/3) = 0, 



it follows from (57), (58), (59), and (60) that Cwy = c^z = Cxy = Cxz = 0, so the 
triangle has no rotation relative to the planes wy,wz,xy, and xz. Since there 
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is no rotation relative to the plane yz either, i.e. Cyz = 0, we have = 0, so 



yz- zy -- 
that 1^ 

at - 



0, and if we assume that z does not take zero values, we can conclude 
= 0, so y{t) = 7z(t), where 7 is a constant. If we denote e = 1 — 7^, 



system ( 63 ) reduces to the family of first-order systems 



(64) 




2m(5-9€'^z*) 



22-1)2 (3e22 + l)2 



Notice that since a point {w, x, y, z) on satisfies the equation w'^+x^+y'^ — z"^ = 
— 1 and z > 1, we necessarily have that e > 1. 

Standard results of the theory of ordinary differential equations can now be 



applied to system ( 64 ) to prove the existence and uniqueness of analytic negative 



elliptic Lagrangian rotopulsators, for admissible initial conditions. To show that, 
except for a negligible set, they are all rotopulsators, we identify the relative 



equilibria, which are fixed points of system (64). 



One fixed point of the vector field in (64) is obviously {z,u) = (0,0), but it lies 



outside the domain z > 1. The other fixed points, if any, must be of the form 
{z, 0), where the positive values of z are given by the roots of the polynomial 

Q{z) = 27{h^ - W)e'z^ - 18{h^ - Ibm^^z^ - Sh^ez^ - - 75m\ 

By Descartes's rule of signs, we must distinguish between two cases: 

(i) \h\ < y/TEm'^, when Q has no positive roots at all; 

(ii) \h\ > vTSm^, when Q has exactly one positive root. 

Case (i) leads to no fixed points for system (64). Case (ii) always leads to one 
fixed point since the unique positive root, Zq := ZQ^m, h, e), is also larger than 1. 

Then we can associate to Q the polynomial 



ez 



To show this fact, let ^ 

g(0 = 27{h' - 9m')e 



Since e > 1, to prove that 2:0 > 1 it is enough show that ^0 '■= e-^o > 1- For 
this, notice that Q{1) = — 48m^ < 0. Since lim^_>.+oo Q(0 = +00, Q has a single 
positive root that is also larger than 1. This remark completes the proof. □ 



11. Negative hyperbolic rotopulsators 

In this section we analyze the solutions given in Definition |4j We first prove a 
criterion for finding such solutions, then provide the conservation laws, and finally 
discuss a particular class of examples, namely the negative hyperbolic Eulerian 
rotopulsators of the 3-body problem in M.^. 

11.1. Criterion for negative hyperbolic rotopulsators. The following result 
provides necessary and sufficient conditions for the existence of negative hyper- 
bolic rotopulsators in H^. 



Rotopulsators of the curved A'^-body problem 29 



Criterion 4. A solution candidate of the form ([8j) zs a negative hyperbolic ro- 
topulsating orbit for system Q if and only if c^x = 0, 

Ej=i"^iK + 1) 

where a is a constant, and the variables Wi, Xi, i = 1,2, . . . , N, satisfy the system 
of 2N second-order differential equations 



(66) 



where 



{wiXi - XiWiY + wl + xf o^{wf + x1 + 1) 



(67) Hi := ^ ' ' ' '\ • -^ • ^ 

i = 1,2, . . . , N , and, for any i,j G {1, 2, ... , A^}, z/jj is given by 

(68) Uij := Qj □ c|j = WiWj + XiXj — {wf + + 1)^ (w^ + + 1)2 cosh(6j — bj). 



Proof. Consider a solution candidate of tlie form ([8j) subject to tlie above initial 
conditions. Then, for any i,j G {1, 2, ... , A^} we have 

Uij := qi □ Qj = WiWj + XiXj - {wj + x] + l)^(w| + + 1)^ cosh(6i - bj), 
and for any 2 = 1, 2, . . . , we find that 

. ^ . iwiXi - XiWiV + wf + xj + {wf + xf + 1)^/3^ 

For all z = 1, 2, . . . , A^, each rj can be expressed in terms of yi and Zi to obtain 

?7i = (Wi +Xi + 1)2, ^i = -^ — T, 

.. _ {wf + xf + l){WiWi + XjXj) + wf + xj + {WjXj - XjWiY 

{wj + xj + l)i 

Substituting a solution of the form ([s]) into system ^ and employing the above 
formulas, we obtain for the equations corresponding to Wi and x, that 

(69) w. = Y: ^ + ^^"^"^ " "t'^t ^ "-^"^ + (^.^+x.^ + l)^^./3^ 



/'7n^ _ 'm.j(xj + z/jjXj) [(w^Xi x^w^) + + xjx^ o , ^2 , 1 x 02 
UUj X, - 2^ ^^^^^ _ + ^2 + ^2 + 1 + l^. + 
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whereas for the equations corresponding to jji and Zi, after some long computations 



that also use (69) and (70), we are led either to identities or to the equations 

N 



(71) 



■A A mjiij sinh(6,- — L 

[Uf- — 1)2 



-, 1 = 1,2,. ..,N. 



For every i = 1,2, . . . , N , we further multiply the ith. equation in (53) by rriiri 
add the resulting N equations, and notice that 



N N 

EE 

i=l j=l 

Thus we obtain the equation 

N 



niimjrjirij sinh(6j — hi 
(4-1)^ 



0. 



i=\ ' ^ i=l ' 



which has the solution 

/3 = 



where a is an integration constant. Consequently, equations (69) and (70) become 
(72) 



N 



rrijiwj + UijWi) [{w^Xi - XiWif + wf + xf\Wi 



(4-1)^ 



wf + xj + 1 



a^iwf + xf + l)wi 



(73) 



Xi 



mj{xj + PijXi) [{wiXi - XiWiY + wf + xf\Xi a^iwf + xf + l)xi 



N 



wf + xf + 1 



i = 1,2, . . . ,N . Since equations (72) and (73) lead to system (66), and the 
condition c^x = occurs from Definition |4| the criterion is proved. □ 

Remark 16. Notice that the constant Cyz of the angular momentum is 

N N 



C. 



yz 



^ mi{yiZi - ijiZi) = mi{wf + xf + 1) = -a. 



i=l 



1=1 



Rotopulsators of the curved A^'-body problem 31 

11.2. Conservation laws for negative hyperbolic rotopulsating orbits. 

In addition to Criterion [4], we would also like to obtain the conservation laws 
specific to negative elliptic rotopulsating orbits. They follow by straightforward 
computations using the above proof, the integral of energy, and the six integrals 
of the total angular momentum. 

Proposition 7. If system ^ has a solution of the form then the following 
expressions are constant: 
— energy, 



i=l 



^ " • -^2 1 „;,2 I ^21 

2{w^ + xf+l] 



mi[{wiXi - XiWiY + wf + xf\ 



(74) 



— total angular momentum relative to the wx-plane, 

N 

(75) c^a; = ^ mi{wiXi - WiXi); 

i=l 

— total angular momentum relative to the wy-plane, 

rnAxiiwiXi — XiWi) — Wi\ . , ,^ 

^ ; ; , — - smh(/3 + k) 

(76) ' \ ' 

Ej=i"^jK + a;j + l) ^ 

— total angular momentum relative to the wz-plane, 

= y ^^^^^^ + 

^ ' N 

+ =^]v o — 7^ ^i'^i (^*^ + a;^ + 1) 5 sinh(/3 + 6^) ; 

Ei=i"^iK + a;- + 1) tl^ 

— tota/ angular momentum relative to the xy-plane, 

N 

(78) 



iV 



^mjXi(wf + + 1)2 cosh(/3 + 6j); 



32 Florin Diacu and Shima Kordlou 

— total angular momentum relative to the xz-plane, 

niAwiiwiXi — WiXi) — Xi] , ,^ 
^ , , ;/i cosh /3 + k 



(79) 



N 

a 

4 



o — Try^'^i^'^^^i + + ^^^^^^^ + 



11.3. Negative hyperbolic Eulerian rotopulsators. We further provide a 
class of specific negative hyperbolic rotopulsators of the curved 3-body problem, 
namely Eulerian orbits in H^. These systems rotate relative to the yz-plane, but 
have no rotations relative to other base planes. 

Consider three equal masses, mi = m2 = =: m, and a candidate solution of 
the form 

(80) q = (qi,q2,q3), qi = iwi,Xi,yi,Zi), z = 1,2,3, 

wi = 0, xi = 0, yi = sinh/3(t), zi = cosh/3(t), 

Wi = w(t), Xi=x(t), yi = ri{t) sinh (3 (t), Zi = ri{t) cosh (3 (t), 

Wi = —w(t), xi = —x(t), yi = ri{t) sinh (3 (t), Zi = rjit) cosh (3 (t), 

with tf ^ + — r]^ = —1. With the help of Criterion |4[ we can now show that, in 
general, these are indeed solutions of system (|4]). 

Proposition 8. Consider the curved 3-body problem in given by system Q. 
Then, except for a negligible set of solutions formed by negative hyperbolic Euler- 
ian relative equilibria, every candidate solution of the form (80) is a negative 
hyperbolic Eulerian rotopulsator. 



Proof. Let us consider a candidate solution of the form (80 ). Then, using Criterion 
|4} straightforward computations show that 

1^12 = 1^13 = {w^ + x^)^ 1/23 = -2(^2 + x^), /3 



m{2w^ + 2x2 + 3) ' 



and the equations of motion reduce to the system 

w = J{w, X, w, x)w 



(81) 
where 

J{w, X, w, x) 



X = J{w, X, til, x)x, 

m(w2 + x2)5 777, 



(w2 + x2 - l)i {2w^ + 2x2 + 1)2 (2ti;2 + 2x2 - l)i 



(WX — XW)2 + ^2 + a;2 _^ ^2 _^ -^^ 



w2 + x2 + l m2(2w2 + 2x2 + 3)2' 
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Notice that the energy integral (74) takes the form 
m[{wx — xwY + w'^ + x"^] 



w"^ + x^ + 1 



+ 



2m(2w2 + 2x2 + 3) 



+ 



2rrr{w^ + x^ 



(w2 + - 1)5 (2^2 + 2x2 + l)i(2«;2 + 2x2 - 1)5 

which imphes that J can be written as 

h m[2(M;2 + x2)(2w;2 + 2x2-1) + 11 

J(uj x) = I ^-^ — - - 

m (2«;2 + 2x2 + l)5(2«;2 + 2x2- 1)1 



m(w2 + x2)5(2^i;2 + 2x2 - 3) 

(^2 + x2 — 1)1 



2m2(2w2 + 2x2 + 3)2' 



From (81), we can conclude that ivx = wx, which implies that 



wx — wx 



k (constant). 



But since from (75) we have that 2m{wx — wx) = c^jx, it follows that k 



From (76), (77), (78), and (79), we can conclude that c 



^10 J/ 



2m ■ 



so the orbit has no rotation relative to the planes wy, wz, xy, and xz. Since there 
is no rotation relative to the plane wx either, i.e. c^ix = 0, we have A; = 0, so 
wx — wx = 0, and if we assume that x does not take zero values, we can conclude 
that 4:^ = 0, so w{t) = ^x{t), where ,^ is a constant. If we denote ( = ^'^ + 1 > 1, 



dt x 



equations (81) reduce to the family of first-order systems 



(821 



u 



h_ I m(4C^x^-2Cx^+l) 
m ^ (2Ca;2+l)l/2(2Cx2-l)3/2 



m.(CJ^)l/^(2Ca:^-3) 
(Cx2-l)3/2 



2m2(2C3;2+3)2 



X. 



Standard results of the theory of differential equations can now be applied to 



system (82) to prove the existence and uniqueness of analytic negative hyperbolic 



Eulerian rotopulsators, for admissible initial conditions. The relative equilibria 



that occur from fixed points of system (82) are given by polynomial equations, so 
they form a negligible set. This remark completes the proof. □ 



12. Negative elliptic-hyperbolic rotopulsators 

In this section we analyze the solutions given in Definition [5j We first prove a 
criterion for finding such solutions, then provide the conservation laws, and finally 
discuss a particular class of examples, namely the negative elliptic-hyperbolic 
Eulerian rotopulsators of the 3-body problem in H^. 
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12.1. Criterion for negative elliptic-hyperbolic rotopulsators. The follow- 
ing result provides necessary and sufficient conditions for the existence of negative 
elliptic- hyperbolic rotopulsators in H^. 

Criterion 5. A solution candidate of the form ^ is a negative elliptic-hyperbolic 
rotopulsator for system ^ if and only if 

(83) d = , P = Tf , 

2^ rriirf M + rriirf 

with di,d2 nonzero constants, and the variables r^, i = 1,...,N, satisfy the N 
second-order differential equations 

^2 



+ 



l + rf 



(84) mj[rj{l + rf) cos(aj — Oj) — rj(l + rf)2(l + r|)2 cosh(6j — hj)] 

where, for any i,j G {1,2,..., A^} with i ^ j , we denoted 

6ij := □ qj = TiTj cos(aj - aj) - (1 + rl)^(l + r|)2 cosh(6i - hj). 

Proof. Consider a candidate solution of the form (|9| for system (|4]). By expressing 
each pi in terms of r^, z = 1, 2, . . . , A^, we obtain that 

r/, = (l + r,)2, 71i= V^ = 77—^73 . 

(l + rf) 2 (l+r/)2 

6ij := Qi □ qj = TiTj cos(aj - aj) - (1 + r,^)2(l + r|)5 cosh(6i - bj), 

q, □ q, = + r^^d^ - - (1 + r^/?'- 

l + rf 

Substituting a solution candidate of the form ^ into system Q, and using 
the above formulas, we obtain for the equations corresponding to Wi and Xi the 
equations 

. 2 

r, = r,(l + r2)(d2_/32)^ 



(85) _|_ j.2^ cos(aj — Oj) — rj(l + ) 2 (1 + r|) 2 cosh(6j — bj 

(86) r,ci + 2r,d = V^^— z = 1, 2, . . . , iV, 
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whereas for the equations corresponding to jji, Zi, we find equations (85) again as 
well as the equations 



N 



^7) 



+ 2fi/3 



i^l - 1) 



1,2, 



iV. 



We can solve equations (86) the same way we solved equations (15) and obtain 

di 



a 



N 2 ' 

i=l ^i^i 



where di is an integration constant. To solve equations (87), we proceed similarly, 
with the only change that for each i = 1,2, . . . , N, the corresponding equation 



gets multiplied by mj(l + rf)2 instead of m^rj, to obtain after addition that 



/3 



do 



M + Ef=i m 



.2' 



where M = X]i=i ^"^^ '^2 is an integration constant. Then equations (85) take 
the form ([84j), a remark that completes the proof. □ 



a 



N 



i=l 



di. 



Remark 17. Notice that the constants c^^^ and Cyz of the total angular momentum 
axe 

N N 
i=l 

N 

Cyz = ^ miiViZi - Ziyi) = "/^ X] ^"-^^ + ^"^^^ ^ 

i=l i=l 

respectively, so di and ^2 describe the rotation of the particle system relative to 
the wx-plane and ?/2;-plane, respectively. 



Remark 18. From (|83|), we can conclude that a and /5 are not independent of 
each 

(88) 



each other, but connected by the relationship 

d2 di 

13 a 

The following observation refers to all the five criteria proved in this paper. 



M. 



Remark 19. There are 2 first-order and 2N second-order differential equations in 
Criteria [2] and |5j but only 1 first-order and N second-order differential equations 
in Criteria [1} [3} and|4j This fact, however, leads to no contradiction because the 
position and velocity of each body rrii in §^ or is given by 6 scalar variables, 
namely {a, P,ri,a, P,ri) for Criteria [2] and [sj {P,Wi,Xi, /3,Wi,Xi) for Criterion |4, 
and {a,yi, Zi,a,yi, Zi) for Criteria [l] and |3| in agreement with the statements o" 
these results. 
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12.2. Conservation laws for negative elliptic-hyperbolic rotopulsators. 

In addition to Criterion |5| we would also like to obtain the conservation laws 
specific to negative elliptic rotopulsating orbits. They follow by straightforward 
computations using the above proof, the integral of energy, and the six integrals 
of the total angular momentum. 

Proposition 9. If system ^ has a solution of the form then the following 
expressions are constant: 
— energy, 

^ rriirf d\ d\ mirrijSij 



^(]\ h = S~^ ' ' I ^ I ^ k 



^2(l + rD 2Ef=im,r| 2(M + Ef=, rn,r|) ,<^^<^ - 1) 
total angular momentum relative to the wx-plane, 



1 1 

2 



N 

Cwy = — 2^ ~~L cos(a + ttj) sinh(/3 + hi 



^ ii 

(90) ^ 

+ mjrj(l + r-j^) 2 [d sin(Q; + Oj) sinh(/3 + 6j) + /3 cos(a + Oj) cosh(/3 + hi)]] 

i=l 

— total angular momentum relative to the wz-plane, 



N 

miTi 



Cu,2 = - cos(a + ai) cosh(/3 + h) 

(91) ^ 

+ ^mjrj(l + rlY\dL sin(Q; + a,) cosh(/3 + 6i) + /3 cos(a + a^) sinh(/3 + 



j=i 

total angular momentum relative to the xy-plane, 

N 



Cxy 



y] sin(a + ai) sinh(/3 + hi) 

r^. 1 + r2 2 



(92) ^ ^ 
+ mjrj (1 + Tj^) 2 [/3 sin(a + Oj) cosh(/3 + 6i) — a cos(a + Oj) sinh(/3 + ; 

— tota/ angular momentum relative to the xz-plane, 

N 

TTLiVi 

(93) „ -'"^ 



sin(a + ctj) cosh(/3 + hi) 

+ "iniTiiX + r^^) 2 [/3 sin(a + Oj) sinh(/3 + 6j) — d cos(a + Oj) cosh(/3 + 
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12.3. Negative elliptic-hyperbolic Eulerian rotopulsators. We further pro- 
vide a class of specific negative elliptic-hyperbolic rotopulsators of the curved 3- 
body problem in M.^, namely Eulerian orbits. These systems rotate relative to the 
wx and yz planes, but have no rotation with respect to the other base planes. 

Consider three equal masses, mi = m2 = := m, and a candidate solution of 
the form 



(94) 



q = (qi,q2,q3), qi = {wi,Xi,yi,Zi), i = 1,2,3, 

Wi = 0, xi = 0, yi = sinh/3(t), zi(t) = cosh/3(t). 



W2 = r{t) cos a(t), X2 = r{t) sin a{t), y2 = rjit) sinh (3 (t), Z2it) = rjit) cosh (3 (t), 
W3 = —r(t) cos a{t), X3 = —rit) sma{t), y^ = rjit) sinh. (3 (t), Z3 = rjit) cosh (3 (t), 

with a and /3 nonconstant functions and r"^ — yf = —1. We can now prove a result 
which shows that, in general, these orbits exist. 

Proposition 10. Consider the curved 3-body problem in given by system 
Then, except for a negligible set of solutions formed by the negative elliptic- 
hyperbolic Eulerian relative equilibria, every candidate solution of the form (94) 



is a negative elliptic-hyperbolic Eulerian rotopulsator, which rotates only relative 
to the planes wx and yz. 



Proof. Consider a candidate solution of the form (94) of system (Q. Then the 
variables relevant to Criterion [5] take the form 



(95) 



>12 



a 



621 = S 



13 



^31 



-(1 + r 



2x1/2 



^23 = 



32 



di 



2mr2 



/3 



m(3 + 2r2) 



with di,d2 7^0 constants. 



In system (84), the equation corresponding to fi is an identity, whereas the equa- 



tions corresponding to r2 and rs are identical, given by 



m(5 + 4r2 



X _l_ ^2 4r2(l -)- ^2^1/2 • 
Therefore, the problem reduces to the family of first-order systems 



(96) 



= P 
p = r(l + 



4 



4 



Am?r^ m2(3+2r2)2 



+ 



rp^ m(5+4r^) 

l+r2 4r.2(i+r2)i/2 ■ 



Standard results of the theory of differential equations can now be apphed to 



system (96) to prove the existence and uniqueness of analytic negative hyperbolic 



Eulerian rotopulsators, for admissible initial conditions. The fixed points of sys- 



tem (96) are obtained from equations that reduce to polynomials, so the relative 
equilibria form a negligible set. This remark completes the proof. □ 
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13. Qualitative behaviour of the negative rotopulsators 

In this section we will prove a result that describes the qualitative behaviour of 
rotopulsators in H^. For this purpose, we first introduce a geometric topological 
object that plays in the role the Clifford torus plays for us in S^, and briefly 
present some of its properties. 

13.1. Hyperbolic cylinders. The 2-dimensional manifold defined by 

(97) Crr, ■■= {iw,x,y,z) eM^'^ \ -7]"^ = -1, < < 27r, ^ G M}, 

where w = r cos 6,x = r sin 6,y = rj sinh ^, and z = 1] cosh ^, with r, p > 0, is called 
a hyperbolic cylinder, and it has constant positive curvature for r and t] fixed. 
But Crrf also lies in M.^ because the coordinates w, x, y, z satisfy the equation 

+ x'^ + y^ - z"^ = -1. 

When r (and, consequently, rj) take take all admissible positive real values, the 
family of hyperbolic cylinders thus defined foliates H^. 

We have previously shown that negative relative equilibria rotate on hyperbolic 
cylinders, [B], [7]. We will next prove that, at every moment in time, a negative 
rotopulsator passes through a different hyperbolic cylinder of any given foliation 

of e^. 

13.2. The qualitative behaviour of negative rotopulsators. We can now 

state and prove now the following result. 

Theorem 2. Consider a negative elliptic, negative hyperbolic, or negative elliptic- 
hyperbolic rotopulsator of the curved N-body problem in . Then, for the foliation 
{Crr])r,ri>o of M.^ givcn by hyperbolic cylinders, the trajectory of each body is not 
contained in a single hyperbolic cylinder, but intersects a continuum of hyper- 
bolic cylinders. For negative elliptic, negative hyperbolic, and negative elliptic- 
hyperbolic relative equilibria, however, each body moves on some hyperbolic cylin- 
der of the foliation for all time. 

Proof. The qualitative behaviour of negative relative equilibria in H'^, including 
more specific properties, was previously described and proved to take place, [6], 
[7j, so we will further focus on negative rotopulsators. So consider the foliation 



[Crr])r,T]>Q of with hypcrbolic cylinders, where Crri is defined in (97). The 
negative elliptic-hyperbolic rotopulsators are given by ([s]), with r and i] noncon- 
stant functions satisfying the relation r"^ — = —1. Each hyperbolic cylinder 
in the foliation, however, corresponds to a constant r and a constant rj that also 
satisfy the relation r"^ — rj^ = —1. Consequently, at every moment in time, each 
particle passes through a different hyperbolic cylinder, thus travelling through a 
continuum of such objects during its motion. 
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In the case of negative elliptic rotopulsators, the rotation takes place only with 
respect to the ifx-plane, so we can treat this situation as a particular form of 
negative elliptic-hyperbolic rotopulsators, for which the result was proved above. 
Similarly, the negative hyperbolic case can be viewed as a particular for negative 
elliptic-hyperbolic rotopulsator, for which the rotation takes place only relative to 
the yz-plsnae. This observation completes the proof. □ 

14. Final remarks 

Criteria 1, 2, 3, 4, and 5 open the possibility to study large classes of rotopul- 
sators for specific values of > 3. Since the examples we provided in this paper 
are restricted to = 3 and regular configurations, it would be interesting to 
know whether there exist rotopulsators given, for instance, by irregular polygons 
for equal or nonequal masses. Other geometric shapes could also be investigated. 

We derived the equations of motion for the positive elliptic, positive elliptic- 
elliptic, and negative elliptic Lagrangian rotopulsators and for the negative hyper- 
bolic and negative elliptic-hyperbolic Eulerian rotopulsators in the 3-dimensional 
case, but used them mainly to obtain existence and uniqueness results. It would 
be interesting to study these equations in detail and provide a complete classi- 
fication of these motions, as it has been done in [9] for the 2-dimensional case. 
Since the equations for these orbits in and are more complicated than the 
ones in and H^, each system would require and extensive study. The stability 
of these orbits, using the tools developed in [H] , is another topic that merits close 
attention. 

Finally, the nonstandard relative equilibria, presented for the first time in this 
paper, deserve special attention themselves. It would be interesting to identify 
new classes of such solutions and understand their behaviour. 

Acknowledgements. Florin Diacu is indebted to Carles Simo and James Mon- 
taldi for some enlightening discussions. The authors also acknowledge the partial 
support of an NSERC Discovery Grant of Canada (Florin Diacu) and a University 
of Victoria Fellowship (Shima Kordlou). 

References 

[1] J. Bertrand, Theoreme relatif au mouvement d'un point attire vers un centre fixe, C. R. 

Acad. Sci. 77 (1873), 849-853. 
[2] W. Bolyai and J. Bolyai, Geometrische Untersuchungen, Hrsg. P. Stackel, Teubner Verlag, 

Leipzig-Berlin, 1913. 

[3] J.F. Carifiena, M.F. Ranada, and M. Santander, Central potentials on spaces of constant 
curvature: The Kepler problem on the two-dimensional sphere and the hyperbolic plane 
H2, J. Math. Phys. 46 (2005), 052702. 

[4] F. Diacu, On the singularities of the curved A'^-body problem. Trans. Amer. Math. Soc. 
363, 4 (2011), 2249-2264. 



40 Florin Diacu and Shima Kordlou 

[5] F. Diacu, Polygonal homographic orbits of the curved 3-body problem, Trans. Amer. Math. 

Soc. 364 (2012), 2783-2802. 
[6] F. Diacu, Relative equilibria in the 3-dimensional curved 7i-body problem. Memoirs Amer. 

Math. Soc. (to appear). 
[7] F. Diacu, Relative Equilibria of the Curved N-Body Problem, Atlantis Press, 2012. 
[8] F. Diacu, R. Martinez, E. Perez-Chavela, and C. Simo, On the stability of tetrahedral 

relative equilibria in the positively curved 4-body problem, |arXiv:1204.572 9 (32 pages). 
[9] F. Diacu and E. Perez-Chavela, Homographic solutions of the curved 3-body problem, J. 

Differential Equations 250 (2011), 340-366. 
[10] F. Diacu, E. Perez-Chavela, and J. Guadalupe Reyes Victoria, An intrinsic approach in 

the curved A'^-body problem. The negative curvature case, J. Differential Equations 252 

(2012), 4529-4562. 

[11] F. Diacu, E. Perez-Chavela, and M. Santoprete, The A^-body problem in spaces of con- 
stant curvature. Part I: Relative equilibria, J. Nonlinear Sci. 22, 2 (2012), 247-266, DOI: 
10.1007/s00332-011-9116-z. 

[12] F. Diacu, E. Perez-Chavela, and M. Santoprete, The A^-body problem in spaces of con- 
stant curvature. Part II: Singularities, J. Nonlinear Sci. 22, 2 (2012), 267-275, DOI: 
10.1007/s00332-011-9117-y. 

[13] W. Killing, Die Mechanik in den nichteuklidischen Raumformen, J. Reine Angew. Math. 
98 (1885), 1-48. 

[14] V. V. Kozlov and A. O. Harin, Kepler's problem in constant curvature spaces, Celestial 
Meek. Dynam. Astronom 54 (1992), 393-399. 

[15] H. Liebmann, Die Kegelschnitte und die Planetenbewegung im nichteuklidischen Raum, 
Berichte Kdnigl. Sdchsischen Gesell. Wiss., Math. Phys. Klasse 54 (1902), 393-423. 

[16] H. Liebmann, Uber die Zentralbewegung in der nichteuklidische Geometric, Berichte 
Kdnigl. Sdchsischen Cesell. Wiss., Math. Phys. Klasse 55 (1903), 146-153. 

[17] N. I. Lobachevsky, The new foundations of geometry with full theory of parallels [in Rus- 
sian], 1835-1838, In Collected Works, V. 2, GITTL, Moscow, 1949, p. 159. 

[18] J.E. Marsden, Lectures in Mechanics, Cambridge University Press, 1992. 

[19] J. Montaldi, Relative equilibria and conserved quantities in symmetric Hamiltonian sys- 
tems, in Peyresq Lectures on Nonlinear Phenomena (Peyresq 1998/1999), R. Kaiser and J. 
Montaldi, eds., pp. 239-280, World Scientific, 2000. 

[20] E. Pcrez-Chavela and J.G. Reyes Victoria, An intrinsic approach in the curved A^-body 
problem. The positive curvature case. Trans. Amer. Math. Soc. 364, 7 (2012), 3805-3827. 

[21] E. Schering, Die Schwerkraft im Gaussischen Raume, Nachr. Kdnigl. Cesell. Wiss. 
Cdttingen 13 July, 15 (1870), 311-321. 

[22] A.V. Shchepetilov, Nonintegrability of the two-body problem in constant curvature spaces, 
J. Phys. A: Math. Cen. V. 39 (2006), 5787-5806; corrected version at math.DS/0601382[ 

[23] A. Wintner, The Analytical Foundations of Celestial Mechanics, Princeton University Press, 
1947. 



